Decoherence-free preparation of Dicke states of trapped ions by collective stimulated 

Raman adiabatic passage 
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We propose a simple technique for the generation of arbitrary-sized Dicke states in a chain of 
trapped ions. The method uses global addressing of the entire chain by two pairs of delayed but 
partially overlapping laser pulses to engineer a collective adiabatic passage along a multi-ion dark 
state. Our technique, which is a many-particle generalization of stimulated Raman adiabatic passage 
(STIRAP), is decoherence-free with respect to spontaneous emission and robust against moderate 
fluctuations in the experimental parameters. Furthermore, because the process is very rapid, the 
effects of heating are almost negligible under realistic experimental conditions. We predict that 
the overall fidelity of synthesis of a Dicke state involving ten ions sharing two excitations should 
approach 98% with currently achievable experimental parameters. 

PACS numbers: 03.67.Mn; 03.67.Lx; 03.67.-a 



I. INTRODUCTION 

A clear understanding of many-body entanglement is 
centrally important for the fundamental description of 
microscopic systems. Within the framework of quan- 
tum information science, entanglement may be viewed 
as a resource for the processing of information in ways 
not permitted by classical logic This has stimulated 
an intense research effort aimed at studying the proper- 
ties of multipartite entanglement, and in turn, experience 
gained from such studies is providing fruitful insights into 
the complex behaviour of condensed matter systems Q • 
However, the characterisation and quantification of mul- 
tipartite entanglement remains an open problem; a full 
classification has only been achieved for small numbers of 
qubits 0, H H, S, 0] and the extension to larger systems 
represents a formidable challenge. 

The experimental generation and investigation of mul- 
tipartite entangled states is therefore a highly desirable 
tool for the future characterisation of many-body entan- 
glement. In this regard, an interesting class of such states 
are the Dicke states, \W^), in which to excitations are 
distributed evenly amongst N parties in an equal coher- 
ent superposition [8|, |9J]: 



(1) 



Here, the symmetrization operator, S creates an 
(un-normalised) equal superposition of all dis- 
tinct permutations of the N qubits, which number 
= N\/[m\(N — to)!]. Dicke states exhibit genuine 
multi-body entanglement [13, El and are also robust 
against particle loss and measurement [HI, E3] • Further- 
more, the entanglement contained in such states cannot 
be destroyed by local operations performed on any 
constituent particle making them interesting from 
the point of view of quantum communication between 




many participants. The notation |Wj^) emphasises 
that the Dicke states represent a generalisation to to 
excitations of the more ubiquitous W-states, in which a 
single excitation is shared evenly between N sites. We 
note that while W-states of up to ei ght ions have been 
created experimentally in ion traps [15|, the extension 
to higher numbers of excitations represents a major 
challenge, both from a theoretical and experimental 
perspective. This is because the dimension of the 
relevant Hilbert space, H., grows very rapidly with the 
number of excitations involved. 

Theoretical proposals exist for the generation of Dicke 
states in a number of physical systems , includin g en sem- 
bles of neutral atoms [3, Qj], [H[ , trapped ions Il9l l2p| , 
quantum dots [2l[ and using linear optics [3, llc| . Ion 
traps are perhaps the best suited of these for the ma- 
nipulation of entanglement in matter systems since they 
offer an unparalleled level of experimental control; the 
state of individual qubits can be initialised, manipulated 
and read out with > 99% precision, and information can 
be stored without loss for times many orders of magni- 
tude longer than typical gate operation times. By us- 
ing carefully chosen laser pulses, the ions' internal and 
motional degrees of freedom can be coupled, and allow- 
ing different ions to interact with a common vibrational 
mode (the 'bus-mode') then creates an effective interac- 
tion between the ions through which entanglement can 
be generated (22j. In this manner, ent ang ling gates have 
been performed with high fidelity [23 . |24| , and complex 
multipartite entangled states have been created [l5l [25l| . 

The generation of large coherent superpositions in con- 
temporary ion-trap quantum information science is lim- 
ited by two major sources of decoherence. These are: 
(i) spontaneous emission as the ions decay from high- 
lying energy states, occupied transiently during internal 
transitions; (ii) motional heating caused by coupling be- 
tween the ions and fluctuating patch potentials on the 
surface of the electrodes used to trap the ions; this has 
the effect of adding or removing vibrational quanta from 
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the system at random. 

Various techniques have been proposed with the aim of 
limiting the effects of spontaneous emission. Of particu- 
lar note is the STIRAP technique for robust transfer of 
population between internal levels of a single particle (26[ 
which has recently been applied to coherent pop ulation 
transfer in a single ion with over 95% fidelity [27|. Here, 
a dark state of the system is utilised, in which higher- 
energy levels always remain unoccupied. In the adiabatic 
limit (whereby the system remains in the dark state at all 
times), the technique is therefore decoherence-free with 
respect to spontaneous emission. Recently, it has been 
suggested to use variants of the traditional STIRAP pro- 
cess in order to implement two-ion quantum gates (28[ 
and to gen erate arbitrary entangled states of a two-spin 
system |29j . 

Regarding the second source of decoherence mentioned 
above, motional heating of the ion chain by its environ- 
ment appears to be unavoidable in traditional multiple- 
step techniques, and the degradation to the state that 
this causes grows with the time taken for the state prepa- 
ration. Therefore the (traditional) approach of using 
many sequential steps, each of which addresses an in- 
dividual ion, appears unsuitable for the creation of large 
entangled states; the Hilbert space dimension and hence 
the number of steps required by an individual address- 
ing approach increases extremely rapidly with both the 
number of ions and the number of quanta involved. 

A promising route to limiting the number of inter- 
action steps is to address all of the ions simultane- 
ously using common laser pulses (a technique known as 
global addressing). In this way, the time required, and 
hence the effects of heating, can be dramatically reduced. 
Several global addressing schemes have been proposed 
flij , [30l l3ll . |32| and experimentally implemented (25| for 
the generation of entanglement in trapped ion systems. 
An alternative route to limiting the effects of heating is 
to look for schemes which do not require the number of 
motional quanta to be controlled precisely [25|, [3(3, [H| • 
In this situation, modest vibrational heating can be tol- 
erated and spontaneous emission is then the dominant 
source of decoherence. 

To date, it has not been possible to find an approach 
to entangled state generation which integrates the ad- 
vantages of insensitivity to vibrational heating along- 
side a robustness against spontaneous emission. Below, 
we present a scheme which goes some way towards this 
aim, by being fully immune to spontaneous emission and 
largely insensitive to heating effects. We propose to gen- 
erate Dicke states of arbitrary size in a chain of trapped 
ions using global addressing by two pairs of overlapping, 
counter-intuitively ordered laser pulses. With an appro- 
priate choice of the laser parameters, it is possible to 
navigate through the overall Hilbert space and into the 
target state whilst remaining in the multi-ion dark state 
at all times. Because the technique presented in this ar- 
ticle is very rapid, the destructive effects of heating are 
largely circumvented, while the use of a multi-ion dark 



state means that spontaneous emission can be avoided 
completely in the adiabatic limit. Consequently, the fi- 
delity of the state preparation using our technique can 
be extremely high. Moreover, the implementation is ap- 
pcalingly simple: it requires only two pairs of sufficiently 
intense and suitably delayed laser pulses, rather than so- 
phisticated sequences of great number of pulses of precise 
areas, as in traditional approaches. 

The remainder of this article is organised as follows: In 
section[TTl a Hamiltonian is introduced to describe a chain 
of N ions interacting with a pair of common laser pulses, 
tuned on resonance with the first red motional sideband. 
Two important properties possessed by this Hamiltonian 
are considered in section lll C[ and by exploiting the sym- 
metry of the chosen system, it is shown in subsection III Dl 
that the system possesses a unique multi-ion dark state 
which is immune to spontaneous emission. In section Hill 
we show that, remarkably, Dicke states of arbitrary size 
can be generated robustly and efficiently in two steps, 
by adiabatic transfer along this dark-state. During the 
transfer process, the dark-state evolves through a net- 
work of increasingly complex A-particle entangled states. 
The net phase acquired during the state preparation is 
identically zero. 

Section IIVI discusses the conditions required for adi- 
abatic following of the dark state, since this is impor- 
tant to the success of our technique. Other sources of 
decoherence, such as vibrational heating and parameter 
fluctuations are considered in section [Vl and it is found 
that the method presented in the current paper is robust 
against these types of decoherence also. In section I VII 
we summarise our findings. 



II. THEORETICAL MODEL AND DARK-STATE 
STRUCTURE 

A. Definition of the problem 

We consider a chain of N identical ions trapped in a 
linear array and cooled to their ground state. Each ion 
has three relevant internal levels, as shown in Fig. [T] 
The computational basis states, |0) and |1), are encoded 
in two hyperfine or Zeeman sub-levels in the electronic 
ground state of the ion and typically have very good co- 
herence properties [H, [H[ . Coupling between these two 
internal states is achieved by a two-photon Raman pro- 
cess, via the upper-level, |e) using only a single pair of 
pulses (labelled a and b) , each of which addresses the en- 
tire chain. Overall, the laser- frequencies, oj a and tub, are 
tuned to the first red motional sideband of the centre-of- 
mass mode for the transition |0) <-» |1), while cu a is tuned 
near to the first red sideband for the |0) «-» |e) transition, 
i.e. 

U a = LOOe - A - v, (2a) 
uj b = cju - A, (2b) 



where woe and u>i e are the Bohr frequencies of the tran- 
sitions |0) \e) and |1) <-» |e) respectively, v is the 
trap frequency and A is a constant single-photon detun- 
ing. We note that in contrast to existing Raman-coupled 
schemes in ion traps [361 ] . the technique proposed below, 
does not require adiabatic elimination of the state |e) 
in order to limit spontaneous emission and in fact we 
choose A < v. In this limit [and because ultimately there 
will be no decay from |e)], the sideband structure on the 
|0) |e) transition can be resolved. In the Lamb-Dicke 
limit and after making the rotating-wave approximation, 
the Hamiltonian for this system (with h = 1 throughout 
this paper) is: 
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fl a (t) [a\e)j(0\j exp (iAt - i<j>f) 

+a f |0)j (e\j exp (-iAt + itf) ] 
+fi i (t)[|e) i (l| i eoq>(tAt-i$) 

+\l) j (e\ j exp(-i&t + i<t>j)[ 



(3) 



Above, and a are the creation and annihilation op- 
erators for centre-of-mass phonons and f2 a is an effec- 
tive Rabi frequency, defined by fl a = r] a £l a / 'y/N, which 
describes the coupling between the motional and in- 
ternal states of the ions 37]. Ci a and £1;, are the 
(real-valued) bare Rabi frequencies for each laser pulse. 
The single-ion Lamb-Dicke parameter for pulse a is 
rj a = \J hk% cos 2 9 a /2M v, with k a being the laser wave- 
number and 9 a the angle between the trap axis and beam 
direction for laser beam a, and M the mass of one ion. 
Equation ([3|) applies in the Lamb-Dicke limit, for which 
T] a (n + l) 1 / 2 <C 1, where n is the number of phonons 
in the centre-of-mass vibrational mode. Also, the trap 
frequency v must satisfy 



i/» |fi (t)|,|fi 6 (i)|,A. 



(4) 



Because the ions in a Paul trap are spaced unevenly along 
the trap axis, there will always be some nonzero phases 



Xjk a cos 9 a — 7r/2, 
Xjkb cos 6b — 7r/2, 



(5a) 
(5b) 



where Xj are the equilibrium positions of the ions and <f> a £ 
the laser phases at the trap centre at t = 0. However, by 
performing a time dependent phase-transformation: 



Hj (t) = U f (i)#j (t)U(t) - iU f (t) 



8U{t) 
dt ' 



(6) 



with 



U(t) =exp 



N 



iE( A *-^)i^-H 



J'=l 



x exp 



(7) 




FIG. 1: Raman-coupled ion qubit. The computational basis 
states |0) and |1) are hyperfine or Zeeman sub-levels of the 
ion's ground state and have very good coherence properties. 
By contrast, decay from the excited state \e) is typically a 
significant source of decoherence in modern ion-trap quantum 
information processing experiments [38l. [39|. 



the Hamiltonian ([3]) can be put into the following simpler 
form, in which each laser couples equally to all ions: 



\e)Mi+tf\Q)Mi 



+ ^[|%(e| i + |e) J -(l| J -]+A| e ) i (e| j 



(8) 



For clarity, we shall first treat the case of equal couplings, 
given by ([8]), and choose to postpone discussion of the 
effects of the phase-factors, <j>^' until section [V] where 
several other technical considerations are also addressed. 



B. Modelling decay from the upper state [e) 

While the two computational basis states can have very 
long coherence times [up to several minutes [13, HU], de- 
cay from the upper state |e) often has a significant effect 
on the dynamics and in current experiments can even 
dominate the error budget for entangling gates [H, Hi| . 
Therefore, leaving aside for the moment vibrational heat- 
ing of the centre-of-mass mode [we shall return to this 
issue in section [V] , the overall dissipative dynamics is 
governed by the following master equation: 



dm 

dt 



p(t),H(t) +Cp{t)- 



(9) 



Here, C is a Liouvillian super-operator corresponding to 
spontaneous emission from the state |e) into either the 
computational basis states, or other levels, at rate T. If 
the upper state |e) is occupied, even transiently during 
the state preparation, the overall fidelity will decrease, 
due to the dissipative term, Cp(t). By contrast, if the 
state |e) remains unoccupied at all times, then this dissi- 
pative term has no effect on the dynamics, and equation 
can be replaced by the Schrodinger equation for a 
pure state, \ip(t)). 
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Below, we shall describe a method for creating Dicke 
states using global addressing, for which the evolution 
takes place along a multi-ion dark state. In the adiabatic 
limit, the coupling between this dark state and all other 
states in the Hilbert space is small enough that sponta- 
neous emission from the level |e) can be neglected. We 
therefore expect spontaneous emission to have only a per- 
turbative effect on the dynamics and so choose to work 
in terms of a state vector, \ip(t)), and the Schrodinger 
equation, rather than a density operator and equation 
©. The effects of emission from state |e) can still be 
included, simply by making the following replacement: 



A = A - iT. 



(10) 



The above replacement is computationally simpler than 
solving the full master equation, ©, and is physically 
equivalent to the assumption that all spontaneous emis- 
sion from |e) takes place into ionic levels other than |0) 
and |1). We note that the norm of the state- vector calcu- 
lated using (|10|) is not constant, since after introducing a 
complex detuning, the Hamiltonian is no longer Hermi- 
tian. The effects of spontaneous emission are therefore 
manifest through a reduction in the norm of \ip(t)). By 
making the replacement (jTUJ) we are taking a conserva- 
tive approach in the sense that the fidelity calculated 
using this reduced-norm state vector is always less than 
or equal to the true fidelity calculated using ([9]). 



C. Symmetries of the Hamiltonian and reduced 
state-space 

The above Hamiltonian, ||SJ), possesses two important 
properties. The first is that the total number of excita- 
tions, 

JV 

JV=ata + 2(|l) J .<l| j + |e) J <e| i ) ) (11) 

2=1 

commutes with Hi(t) and is therefore a conserved quan- 
tity. If the system is prepared in a state with a 
specific number of quanta, m, then the subspace of 
states containing m quanta is closed upon evolution 
under the Schrodinger equation. However even when 
the system is restricted to the subspace containing ex- 
actly m quanta, the dimension of this subspace is still 

E™ o T,™=o C m-/ 1 C £™~' i > which S rows rapidly with both 
m and N. 

A brute-force approach to solving the dynamics is 
clearly impractical. Fortunately however, Hi(t) pos- 
sesses a symmetry that permits a further simplification 
of the dynamics - the Hamiltonian ([5]) is invariant un- 
der interchange of any two ions. Providing that the sys- 
tem begins in a symmetric state, it will always evolve 
through symmetric states at later times and with an ap- 
propriate choice of basis, the dynamics is confined to a 
closed C^ +2 -dimensional subspace of 7i [52| . We choose 



the notation \W^{pL, e)) for a symmetric state of m ex- 
citations shared between N ions and the bus mode, fi 
represents the number of phonons, and e is the number 
of ions in state |e). This leaves m — \l — e ions in state 
|1) and therefore N — m + [i ions in state |0). Formally, 
\W^(fx,e)) is defined as follows: 



m— fi— e N—m+iJ, 



(12) 



where S creates an (un-normalised) equal superposition 
of all distinct permutations of the ions' internal states, 
and the normalisation coefficient, A™(/i, e), is given by 



e\(m - fj, - e)\(N - m + fj,)\ 

All but three of the states defined by Eq. (fl~2"]) 
are symmetric entangled states of all N ions. The 
exceptions are the state | W^(m, 0)} = |0i . . .0jv)|m) 
which is an m-phonon Fock state, with all of 
the ions in state |0), \W%(0, 0)) = |li . . . l m )|0) and 
\W™(0,m)) = |ei . . .e m )|0) which are symmetric prod- 
uct states of the ions' internal states, with no excita- 
tions in the bus mode. We note that |W™(0, 0)) and 
\W™(0, m )) omv exist if the number of excitations is the 
same as the number of ions, m, that we choose to ad- 
dress. We note also that in this notation, the Dicke state 
defined in Eq. ([T]) (and with zero vibrational quanta) is 
written \W%)\0) = \W% (0,0)). 

The Hamiltonian ijHJ) consists of five terms, each of 
which is symmetric under the interchange of any two ions 
and conserves the number of quanta, N. These have the 
following effects: 



n„(t) A 

2 = 1 



2=1 



<W)) 

A+(/x,e)|lC(/i-l,e+l)), (14a) 



A- (/i,e) w£(/i+l, e -l)>, (14b) 



2 2=1 



2 2=1 



1>2<4-J \W£(ji,e)) 

= Xi(n,e)\W%(n,e-l)), 
e)i{l\j] \W%(ji,e)) 



A 



= \+(fi,e)\W»(!i,e + l)), 
\e)M] |^(M,e)>=eA|W^( M ,e)>, 



N 
2 = 1 



(14c) 

(14d) 
(14e) 
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Wl( M, £+ i) W^E+l) 

(8+l)A 




W>,e-l) W:(n+l,e-l) 



(e-l)A 



FIG. 2: Each symmetric state given in Eq. (I12|l is coupled to 
at most 4 other states as described by Eqs. (|14|l and Q5) . 



with the coupling matrix elements, A^ b , given by 

n„(t) 
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VM£+l)(iV-m + M), (15a) 
v/(M+l)e(W-m + // + !), (15b) 



(15c) 
(15d) 



The above coupling structure is sketched schematically 
in Fig. [21 By symmetry, it is clear that these couplings 
must satisfy A+(/i, e) = A~ (p, — 1, e + 1) and A^(/x, e) = 
\h(fi,e + 1). Also, from considering the boundaries of 
this level scheme (which only involves a finite number 
of states), we have A+(0, e) = A~(^,0) = Aj~(/x, 0) = 
A^(/i, m — /i) = 0. All of these conditions are satisfied by 
Eqs. CG3). 



D. Deriving the dark-state, \Lpo(t)} 

By definition, a dark state is an eigenstate of the 
Hamiltonian which does not decay by spontaneous emis- 
sion. Clearly, such an eigenstate cannot involve the state 
|e) in any of the ions and therefore in the interaction- 
picture defined by -Hj(t), a dark state has zero eigenvalue, 
i.e. 



Y,H jk (t)c k (t)=0 Vj. 



(16) 



Here Ck(t) are the elements of the dark-state vector, 

M*)> = X>(*)l*>, (i7) 

k 

and the index k runs over all states defined in equa- 
tion (|12p. The matrix elements Hjk(t) are defined by 



E = 3 



E = 2 



£= 1 



e = 





m-phonon 
state, W*(m,0) 



FIG. 3: The levels \W£{n,e)) are arranged horizontally by 
e and vertically by fj,. Each state is only coupled to states 
in neighbouring horizontal manifolds and the coupling always 
links states with a different value of e. As in other figures, the 
red couplings are proportional to fl a (t) and the blue couplings 
are proportional to f2t,(i). 



Hjk(t) = (j\Hi(t)\k) . It is certainly not obvious a priori 
that there is such a state satisfying Q16p . However, we 
shall prove that there is one, by first assuming its exis- 
tence and then calculating explicitly the elements Ck{t). 

To facilitate the derivation of \(po(i)}, we choose to 
group the states \ W^(fi, e)) into horizontal manifolds ac- 
cording to e, i.e. the number of ions in state |e). The 
other ions are found in either state |0) or |1), and all 
distinct permutations of the ions' states with an appro- 
priate number of excitations are to be included. The 
reason that we have chosen to express the level scheme 
in this way is as follows. By definition, a dark state can 
only involve the states in the lowest (e = 0) manifold, 
since all other states can decay by spontaneous emission. 
Therefore only this manifold and the states which are 
directly coupled to this manifold are relevant. We thus 
define the set of relevant states, M, as the set of all 
states of the system with either zero or one ion in the 
state |e); i.e. all states \W^(fi, 0)), (^ = 0, . . . ,m) and 
|W£(M)>,(A* = 0,...,m-l). 

The reason that these are the 'relevant' states is that 
the other manifolds are all (i) unoccupied and (ii) only 
coupled to unoccupied manifolds. Therefore, we auto- 
matically have: 

J2 H jk{t)c k (t) = Q j$M, (18a) 
H jk (t)c k (t) = jeM, (18b) 

k4M 



and so Eq. (| 16[) reduces to the much simpler condition 
H jk (t)c k (t) = jeM. (19) 



keM 



Therefore, if a dark state exists, the couplings to and 
between states that lie outside of M. are irrelevant, and 
can be neglected, as shown in Fig. [4j Hence, if a dark- 
state exists within the relevant subspace, M, then this 
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Couplings irrelevant 



^ All these manifolds are empty for a dark state 



Only these levels are occupied 



Couplings irrelevant 




Therefore only this coupling block is relevant for the dynamics 



FIG. 4: Couplings within the relevant subspace M. The dark 
state is composed solely from states in the lowest (e = 0) 
manifold. 



must also be a dark-state of the overall Hamiltoman given 
inEq. ©. 

Upon consideration of the couplings shown in Fig. 2J 
it becomes clear that the relevant subspace does indeed 
possess a dark state, since the states of M. are arranged 
in the configuration found in chain-STIRAP [26| . In par- 
ticular: 

• There are always an odd number of states in Ai; 
m of these form an upper manifold, with detuning 
A, while (m + 1) states comprise a lower-manifold, 
with zero detuning. 

• Coupling between the levels alternates in a chain 
between couplings proportional to f2f,(t) and cou- 
plings proportional to ft a (t). 

For ease of notation, the abbreviations A QjM = A+(/k,0) 
and Ab.p = \£ (fx, 0) are used from this point onwards, 
since within the relevant subspace, both the value of e 
and also the label ± can be omitted without ambiguity. 
Furthermore, after noting that the only states that are 
ever involved in the dark state are Dicke states, we may 
revert to the simpler notation introduced in Eq. ([T]) and 
write IW^JIm) in place of \W^(fi,e)). The first ket 
denotes the internal Dicke state of N ions sharing m — fx 
quanta and the second ket represents a motional Fock 
state of /x phonons. Written in block-matrix form, the 
dark-state condition (1191) takes the form 



/ Ab, 
A&.o A A Qj i 

Xa.l 






\ o o 









A&, m -l 









Afe.m-l 

A 





( co \ 









Cl 




Cm — 1 









V c m / 



=0, 



(20) 



which is satisfied if 

Cfj,\b,ii + C^+lAa^+i =0 V (i = 0, . 



We are now in a position to write down the dark state of 
©: 



where the amplitudes c^(t) satisfy 



<>>'■> =-^-\\Kk{t) n 

fc=0 l=ii+l 



(22) 



(23) 



and the normalisation factor A is given by: 



m [ /x— 1 m 

, E n n ^ 

\ M=0 [k=0 l=n+l 



A = 



We note that in Eqs. (|23|) and (|24|) . the convention 



(24) 



k—r 



(25) 



has been used. The state \<fo(t)) is an eigenstate of 
the total Hamiltonian, (JSj) with eigenvalue zero. More- 
over, it can immediately be seen that the dark subspace 
of H(t) is non-degenerate. Indeed, the recurrence re- 
lation (|2Tj) together with the normalisation condition, 
(ip (t)\(pf)(t)) — 1, defines \ipo(t)) uniquely, up to an arbi- 
trary phase, so there is no other eigenstate of H(t) with 
zero eigenvalue. 



III. GENERATION OF DICKE STATES 

As shown above, the composition of the dark state 
depends upon the ratio of the Rabi frequencies Q a (t) 
and Qb(t). It is crucial for the technique proposed here, 
\ipo(t)) can be transformed smoothly between a product 
state and an entangled Dicke state, simply by a controlled 
manipulation of the Rabi frequencies. Therefore, by en- 
forcing adiabatic evolution, the entire state preparation 
can be performed without the system ever leaving the 
dark state. This process is detailed in the following three 
subsections 



A. Properties of the dark-state 

The properties of the state \ipo) that we will use are: 

|fia| > Pb\ — > Ic^+il < M O = 0,...,m- 1), 
|fia| < Pb\ — > Ic^+il > Ic^l O = 0,...,m- 1), 

and in particular: 

|W^)|0), 



1. (21) 



1™, !>()>} 
|n»|/|n o |->-0 

lim {\tfo)} 



|0...0)|m) 



(27a) 
(27b) 
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Equations (|27a|) and (|27b[) are exactly the properties 
of the dark state used to perform complete population 
transfer in STIRAP and chain-STIRAP, and which have 
been studied extensively [26|, |4(| [H| • Here, however, the 
STIRAP approach is extended to create entanglement 
in a multi-particle system, rather than to transfer popu- 
lation between different levels corresponding to a single 
particle. Drawing on the experience gained from previ- 
ous studies in other STIRAP-type systems, we choose to 
use overlapping Gaussian pulse shapes, as these result 
in a favourable convergence to the adiabatic limit with 
increasing pulse area 



Q (t) = O exp 



Qi,(t) = £!q exp 



-{t + r? 

T 2 



(28a) 
(28b) 



We define the terms forward STIRAP and reverse 
STIRAP as follows: 

• Forward STIRAP (r > 0): When pulse a precedes 
pulse b, we have: 



lim i 

t— >-oo I \l„ 



lim <^ \ = 0, (29) 

and so an initial state |Wj^)|0) will be completely 
transferred into the state |0 . . . 0) |m). In the termi- 
nology of single-particle STIRAP, fl a is the Stokes 
pulse and is the pump pulse. 

• Reverse STIRAP (r < 0): When pulse b precedes 
pulse a, we have: 



lim <^>=n 

t— >— OO I Sif, 



and \(po(t)) evolves from |0...0)|m) before the 
pulses, to |VFm)|0) afterwards. Here the roles of 
the two pulses are reversed: fit, is now the Stokes 
pulse whereas il a is now the pump pulse. 

We note that both forward and reverse STIRAP utilise 
a counter-intuitive pulse ordering; the difference between 
them is simply the direction of the adiabatic transfer. 



B. State preparation by adiabatic passage 

Dicke states of any number of ions and quanta can be 
prepared by applying just two pairs of laser pulses, as 
follows: 

1. The chain is cooled to its vibrational ground state, 
|0), with m of the ions initialised in state |1) and 
the others in state |0). 



2. Simultaneously addressing these m ions on the red- 
sideband, but leaving the other N — m ions in their 
initial state, the pulses given by Eq. (|28|) are ap- 
plied, with t > (i.e. forward STIRAP). This has 
the following effect: 



1 1 ... 1 ... 0) |0) — ► |0...00...0)|m). 



(31) 



N- 



N 



If it is technically demanding to address these m 
ions simultaneously without addressing the other 
N—m ions in the chain, this Fock state preparation 
may also be achieved in a sequential manner, by 
performing stage 2 repeatedly on a single ion. 

3. All N of the ions are addressed simultaneously by 
applying the pulses of Eq. (|28|) again, but this time 
using reverse STIRAP (i.e. r < 0). This guides the 
dark state of the system into a Dicke-symmetric 
state, as required 



|0...0)|m) 



vOlo). 



(32) 



As an aside, we note that at the end of stage 2, the system 
resides in a vibrational number state. A detailed discus- 
sion of the generation of motional Fock states by global 
addressing can be found in [42| , wherein a comparison of 
the heating effects and adiabaticity requirements for Fock 
state generation by single- and multiple-ion addressing is 
also given. 



C. Example of Dicke state preparation 

We have tested the technique described above by nu- 
merical solution of the Schrodinger equation for up to 
10 ions and sample results are shown in Fig. [5l Note 
that the two adiabatic passage steps in the Dicke state 
preparation differ only in the number of ions that are 
addressed [m in the first step and N (N > m) in the sec- 
ond]. Therefore, for simplicity, only the second stage is 
plotted, starting in |0...0)|m) and ending in |W,^)|0). 
Figure O shows the example case of N = 5, m — 2, for 
which the total number of coupled many-particle states is 
dimW = 51. These are grouped in three frames accord- 
ing to £ = 0,1,2, and then according to ^ within each 
frame. States with the same e and /i are equivalent up 
to a permutation of the ions and so for clarity only one 
curve is shown for each value of /j, and e. Therefore, the 
population of only six out of the overall 51 basis states is 
plotted. We draw particular attention to the red curve in 
Fig. [5ji, which shows the population of 1 11000) |0) . This 
state, along with 9 other basis states with the positions 
of the ions permuted, comprise the Dicke state, |Wf )|0), 
which explains the final state population of ~ 0.1 for the 
red curve in Fig. [SJ;. 

In the example shown, around 98.5% of the popula- 
tion remains in the dark-subspace (frame(c)) during the 
transfer process. The slight leakage of around 1.5% into 
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-2 -10 12 -2 -1 1 2 

Time (Units of T) Time (Units of T) 

FIG. 5: Dicke state preparation for the case N — 5, m — 
2. Different curves within frames (a), (b) and (c) corre- 
spond to different values of fi. Frame (d) shows the time- 
dependent overlap with the target Dicke state. Parameters 
are fl T = 50, FT = 2, AT = and r = -0.6T. 



decaying states is due to non-adiabatic couplings. The 
conditions for adiabaticity are discussed in the following 
section. 



IV. ADIABATICITY REQUIREMENTS 



A. Definition of the adiabatic basis 



1. Adiabatic states 



stants. However, the explicit time dependence of Eq. 
(j8]) introduces coupling between the different eigenstates, 

\<P(t))- 

^ = -£<^)l|k'(*)><v(*) 

x exp (i [E v (t') - E v ,{t')] dA . (35) 

In the adiabatic limit, the splittings between the various 
E v (t) are large enough and the variation of Hi(t) is slow 
enough that these couplings have a negligible effect, and 
a system prepared in an adiabatic state remains in that 
state without making transitions. The only effect on the 
state is that it acquires a dynamical phase, f.f E v (t') dt' , 
which for the dark-state, \ipo(t)), is identically zero. 

By contrast, if the splitting between the adiabatic en- 
ergies is too small, then the couplings given in (|35[) can 
cause transitions between the eigenstates of Hj(t). The 
terms on the right-hand side of Eq. ([33]) are therefore re- 
ferred to as non-adiabatic couplings. Such couplings are 
harmful for the technique proposed above, since they can 
cause a leakage of population away from the desired state 
and a concomitant reduction in the fidelity. It is there- 
fore very important to know the conditions necessary in 
order to reach the adiabatic limit - these are studied in 
the remainder of this section. 



The term adiabatic basis refers to the time-dependent 
basis formed by the instantaneous eigenstates of the 
Hamiltonian, (j8]). So far in this paper, we have concen- 
trated on only one of the adiabatic states, namely \<po(t)), 
which has eigenvalue zero. Below, it will prove useful to 
consider the other eigenstates, \<p(t)), and eigenvalues, 
E v (t). In the adiabatic basis, the Hamiltonian can be 
written as 

H I {t) = ^E v (tMt)){<p(t)\ (33) 

<p 

and a general state-vector can be expressed as follows: 
m))=Y,^(t)exp E v (t') dt'^j \<p(t)). (34) 

We note that for non-zero T, the energies E v (t) can 
be complex (with negative imaginary component). This 
makes clear physical sense, since the only non-decaying 
eigenstate of H](t) in this situation is the dark state, 
IVo(*)>. 

2. Non-adiabatic couplings 



B. Adiabaticity conditions 

1. Dependence on A and fio 

In order to remain in the state \ipo(t)) throughout the 
entire process, the energy splitting in the eigenspectrum 
of Hj(t) must always be large enough that the effects 
of non-adiabatic couplings between \cpo(t)) and all other 
eigenstates of Hj{t) are negligible. Since we are inter- 
ested in the adiabatic limit, a perturbative treatment of 
the non-adiabatic couplings is valid, and by this reason- 
ing it is clear that the most significant non-adiabatic cou- 
plings are to those states whose instantaneous energies 
are closest to the dark state energy, i.e. zero. We choose 
the label E\ (t) for the eigenvalue of Hj (t) that is closest 
to zero, and \ipi(t)} for the corresponding eigenstate. 

An analytic derivation of the energy E\ appears to be 
highly demanding. Fortunately however, for the purposes 
of estimating the adiabaticity requirements of the above 
process it is sufficient to determine the functional depen- 
dence of E\ on the Rabi frequency, Qq and the detuning, 
A. In order to do this, we note that E\ is a root of the 
characteristic polynomial: 



For a stationary system [with time- independent Hamil- 
tonian], the coefficients, a v (t) in Eq. (f3"4"|) would be con- 



M 2m +i = 0, 



(36) 
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where the determinant Mi m +\ is defined by 



M 2 



—E v 








A n .i — E, x 









Xb,m-1 A — E, 
A 0l 







(37) 



In appendix[A] it is shown that the solutions of Eq. |36|) 
for which E v ^ can be expressed as P(z) — 0, where 
P(z) is an m^-order polynomial in the dimensionless pa- 
rameter z = E V (E V — A)/Qq. The coefficients in P(z) 
depend on m, N and the time, t, but importantly, they 
are independent of both Slo and A. We are interested 
in a particular root of P(z), corresponding to the energy 
Ei, and choose to write this root as z = 7m (i), i- e - 



1 



4Clfa£(t) . (38) 



Now, the larger the value of E\(t) at any given time, 
the smaller the effect of non-adiabatic couplings. There- 
fore, just as in three-state STIRAP [26[ and other chain- 
STIRAP models [Z(|, adiabaticity can be improved ei- 
ther by increasing the pulse amplitude fig or reducing 
the single-photon detuning A. It is interesting to con- 
sider two specific limits: 

• Zero single-photon detuning: In the case A = 0, 
the adiabatic energies come in sign-conjugate pairs, 
and the non-adiabatic couplings are dominated by 
two states with adiabatic energies 



E x {t) = ±Q Q ^Jt). 



(39) 



• Large single-photon detuning: If A 3> Oo, the adi- 
abatic energy is approximately 



E t (t) 



(40) 



which is smaller than in the resonant case. Larger 
values of A require a higher Rabi frequency or a 
slower variation of Hi(t) in order to ensure adia- 
batic following of the state \ipo(t)). 

In order to minimise the various types of decoherence, 
it is in general desirable to use pulses with the shortest 
possible duration. In addition, in order to suppress the 
effects of spontaneous emission in this many-particle STI- 
RAP, the process must also remain adiabatic. Increasing 
the value of f2o permits the use of shorter pulses while re- 
taining a high level of adiabaticity. The maximum value 
of f2o is restricted by the condition that Qq must remain 
significantly less than the trap frequency in order to sat- 
isfy the rotating-wave approximation and to prevent un- 
wanted excitations in off-resonant vibrational sidebands 



[see Eq. ((4]) and Ref. [37|]. A second route to achieving 
shorter interaction steps whilst remaining adiabatic is to 
use resonant pulses, with A = 0. Use of resonant Ra- 
man coupling represents another significant advantage of 
utilising a dark state, since by being much faster than 
traditional (off-resonant) Raman-coupled schemes, the 
process has a natural insensitivity to heating effects. A 
contour plot which demonstrates the dependence of the 
fidelity on Hq and A [in the presence of heavy losses from 
|e)] is shown in Fig. [5] 



2. Dependence on m and N 

In general, 7^ (t) is a complicated function of both the 
number of excitations and the number of ions. A numer- 
ical evaluation shows that the area ^fi /^ 7m(* / ) dt' 
decreases weakly with m and increases in proportion to 
y/~N for large N. However, we note that a characteristic 
feature of adiabatic passage in chainwise-coupled systems 
is that the departure from adiabaticity is determined pre- 
dominantly by the weakest coupling in the chain [40(. 
For the current system, this is always A&o(*) = ^b(t)/2, 
which is independent of both m and N . Indeed, an an- 
alytic expansion of 7„(t) for resonant pulses, large N, 
and for the first few values of m shows that: 



r < 



T > 




ifi (t)\/iv, 



(41a) 
(41b) 



Therefore, for both forward and reverse STIRAP, non- 
adiabatic couplings are at their most potent during the 
tail of the pulse fib (t) . Since the adiabatic energy E\ (t) is 
independent of both m and N in this region, we conclude 
that the adiabaticity condition for Dicke state prepara- 
tion does not depend strongly on either m or N. 



3. Dependence on V 

It is well-known that the STIRAP process is robust to 
moderate losses, since the adiabatic transfer state \yo{t)) 
does not feature the intermediate state, |e). However, for 
very large decay rates, spontaneous emission can deterio- 
rate the fidelity of the population transfer. This deterio- 
ration takes place through two physically distinct routes: 
(i) In the absence of losses, it is possible for population 
to leave the dark state temporarily via the non-adiabatic 
couplings given in Eq. (I35[) but to return to \(po(t)} by 
the end of the process [53j. However, for V 7^ 0, popu- 
lation which leaves \ipo(t)} may be irretrievably lost, (ii) 
The strength of the non-adiabatic couplings is dependent 
on r, and so increasing the rate of spontaneous emission 
can alter the requirements for adiabatic following. In the 
extreme case, where flo/T — > 0, this can lead to quantum 
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FIG. 6: Fidelity of the Dicke state preparation as a function 
of the Rabi frequency, fio, and single-photon detuning, A, for 
m — 2, N = 5, t — — 0.6T and in the presence of heavy losses 
(FT = 2). These are the same parameters as in Fig. [5] The 
large region of high fidelity at large £IqT and small AT show 
that the procedure is robust to significant fluctuations in the 
experimental parameters. 



overdamping, wherein the system becomes frozen in its 
initial state [431 ]. 

A thorough analysis of the effects of losses in three- 
state STIRAP (which corresponds to a single excitation 
in this model) has been performed for the cases of pure 
dephasing [44j|. decay into the computational basis states, 
|0) and |1) |45|, and decay into other states [43|]. When 
the decay from |e) is into the computational basis only, 
the second loss mechanism given above is dominant and 
the transfer efficiency is found to be [HJ 



Transfer efficiency = 1 — exp 



& t 



(42) 



Conversely, if the state |e) relaxes predominantly into 
states outside of the computational basis, then the first 
relaxation channel given above can also play an impor- 
tant role, and the effects of losses are more pronounced. 
Nonetheless, the relevant parameter for adiabaticity is 
still HqT/T. A quantitative assessment of the effects of 
spontaneous emission expected with current experimen- 
tal parameters is given in section IV Al 



C. Dependence on pulse delay 



From the results shown in subsections If If CI and IfV Bl 
it is clear that for a counter-intuitive pulse ordering, the 
dark state \ipo(t)) connects the initial and target states, 
following the pathway described in Sec. IIIII In this case, 
allowing a significant decay rate has almost no effect on 
the dynamics, providing that the evolution remains adi- 
abatic, since \ipo(t)} is a non-decaying state. For an in- 
tuitive pulse ordering, the situation is more complicated; 




-1 41.5 0.5 1 
Delay, r (units of T) 




-1 -0.5 0.5 1 
Delay, T (units of T) 

FIG. 7: Fidelity of the transfer process vs. the time delay, r 
between the pulses. The solid red curves show the idealised 
dynamics in the absence of spontaneous emission, while the 
blue dashed curves include decay from the upper level at rate 
r. Parameters match those in Fig. [S] i.e. m = 2, N = 5, 
n T = 50, FT = 2 and (a) AT = 0; (b) AT = 30. 



it is possible that some population is transferred to the 
target state, via another, non-dark adiabatic state [4l| . 
In this scenario, admitting a finite decay rate, T, should 
greatly reduce the overall fidelity, since the intermediate 
states are then able to decay. 

In order to test the above reasoning, the fidelity is plot- 
ted against pulse delay r in Fig. [7] for forward-STIRAP. 
Figure [7ji shows the case where the single-photon detun- 
ing, A is zero, and Fig. [7p shows an off-resonant case, 
with AT = 30. Positive r corresponds to a counter- 
intuitive pulse ordering, and in this situation, the trans- 
fer process can tolerate significant variations in r, while 
retaining a transfer efficiency close to 100%. Note that 
this robustness is more pronounced in the resonant case. 
Also, allowing a significant decay rate from the state |e) 
has a negligible effect on the dynamics, which explains 
the close overlap between the solid and dashed curves in 
Fig. [7] for positive r. 

For intuitively-ordered pulses, the transfer efficiencies 
can still be high for certain values of the pulse delay in 
the idealised case of T = 0. However, the system then 
evolves though an adiabatic state which has the potential 
to decay, and therefore the fidelity is very low when a 
significant decay rate is introduced. This explains the 
behaviour of the blue dashed curve for negative r. 



1. Adiabatic elimination of the state \e) 

It is interesting to note that the red curve in Fig. [TJd 
is nearly symmetric with respect to the pulse delay. This 
can be understood through an adiabatic elimination of 
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Eliminate upper levels aiiiabalieally 




IO...0>|m> 



H€>|0> 



4 Transition pathway for 

Intuitive pulse ordering 

FIG. 8: For a large single-photon detuning, A, it is possible 
to eliminate the e = 1 manifold adiabatically. In this case, 
the chain-STIRAP model described by Eq. (|20)l reduces to 
an (m + l)-level "bow-tie" energy pattern. The slope of the 
effective detuning, 5(t), of levels in the reduced state space is 
controlled by the pulse delay, r. It is known that in a bow-tie 
crossing, transfer efficiency does not depend on the sign of the 
chirp rate, and this explains the symmetry of the red curve in 
Fig. [7p with respect to the sign of r. Including a finite decay 
rate breaks this symmetry. 



the states in the e = 1 manifold, which is appropriate 
when the single-photon detuning A is large. In this case, 
the chain-STIRAP structure given in Eq. (|2T)]) reduces 
to an (to + l)-level "bow-tie" coupling, as outlined in 
Fig. [5] Different pulse delays in the original picture are 
mapped onto different chirp rates in the bow-tie picture. 
[The Gaussian time-dependence given by Eq. (f2"8")) results 
in an approximate sech-tanh dependence for the bow-tie 
couplings and detunings [43| - i.e. a multi-level Allen- 
Eberly model 14611 .1 Such a bow-tie coupling scheme was 
discussed in [19J, and arises naturally in the prepara- 
tion of Dicke states using two-level ions and a frequency- 
chirped laser pulse. As shown in [l9j . and sketched in Fig. 
[SI an adiabatic pathway connecting the initial and tar- 
get states always exists in the bow-tie model, regardless 
of the sign of the chirp rate. This explains why the red 
curve in Fig [TJd is symmetric under inversion of the sign 
of t. Admitting a finite decay rate breaks this symmetry, 
since, for negative r, the adiabatic state used for the state 
preparation is a decaying state whereas for positive r this 
adiabatic state is non-decaying. The above behaviour is 
closely analogous to that found in traditional STIRAP 
for the case of three levels [4l[ , here observed in a many- 
particle, multi-level scenario. As such, it provides strong 
confirmation that the dynamics do in fact take place via 
a multi-particle dark state, as claimed above. 



V. TECHNICAL CONSIDERATIONS 

Below, we assess the effects of decoherence for cur- 
rently achievable experimental parameters. We note that 
the effects of spontaneous emission can be reduced by us- 
ing longer pulses (as this increases the adiabaticity of the 
process), while vibrational heating is reduced by using 
shorter pulses. Therefore, in order to maximise the over- 
all fidelity, the pulse duration must be chosen in order 
to achieve a balance between these two mechanisms for 
decoherence. For concreteness, all values quoted below 
refer to the specific example of trapped Ca + ions. 



A. Spontaneous emission 

In order to limit the effects of spontaneous emission, 
the system must remain in the dark state, \(po(t)), at all 
times. For all of the ion species favoured in contempo- 
rary quantum information processing experiments using 
Raman-coupled qubits, this is a non-trivial requirement, 
since the decay rate T greatly exceeds the allowed values 
of the Rabi frequencies, fl a (t) and 0&(t). Therefore, the 
adiabaticity of the state preparation is limited by the ef- 
fects of overdamping [47| , as described in section IIVB 31 
Equation (|4"2")l shows that in order to achieve a trans- 
fer efficiency of 1 — x (with x <C 1), the minimum pulse 
timescale is: 



(43) 



For 43 Ca+, the decay rate is V /2-k ~ 22MHz 36], while 
the maximum permitted value of f2rj is limited by the 
dual requirements that the vibrational rotating-wave 
approximation should remain valid [48] and that un- 
wanted excitations are not created in vibrational modes 
other than the centre-of-mass mode (37]. A conserva- 
tive estimate for flo which satisfies both of these re- 
quirements is fio v/lQ WA - Taking a trap-frequency 
of z//2vr ~ 4MHz H, then gives T ~ 80^s in order to 
achieve a transfer efficiency of 99%. 



B. Heating effects 

In order to limit the effects of heating, it is desirable to 
carry out the state preparation as fast as possible, and 
in this sense the fact that our technique requires only 
two interaction steps represents a significant advantage. 
Ultimately, however, Eq. (|4"3"|) sets a lower limit to the 
duration of each interaction step. 

In order to estimate the effects of motional heating 
we assume that the heating rate per ion is ~ 5Hz, as 



recently measured in Ref. [24j . Although this value was 
measured for the isotope 40 Ca + , we note that the same 
value is expected for 43 Ca + (3f| ■ Allowing an overall time 
of ?bicke ~ 6T for the two pairs of pulses, we therefore 
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predict the total number of heating events during the 
state preparation to be on the level of 2.4 x 10 -2 for a 
chain of 10 ions: i.e. an overall fidelity approaching 98%. 

Finally, we note that the effects of heating can be re- 
duced by several orders of magnitude if another vibra- 
tional mode is used, rather than the centre-of-mass mode. 
In this case, individual addressing of the ions is necessary, 
with appropriately tuned couplings, in order to realise the 
equal-coupling Hamiltonian (jSJ) . An auxiliary ion may be 
incorporated in the chain for the purposes of continuous 
sympathetic cooling of the centre-of-mass mode [49| . Be- 
cause the heating rates for higher-order modes are con- 
siderably lower than for the centre-of-mass vibrational 
mode, the effects of motional heating may be greatly 
suppressed in this case. Use of higher-order modes there- 
fore permits an increase in the pulse timescale T which 
helps adiabaticity and hence further reduces the effects 
of spontaneous emission. 



C. Uneven spacing of the ions 

In the main body of this article, we have considered 
a Hamiltonian, Hj(t) [given in Eq. ©], which is sym- 
metric under interchange of any two ions. Now, as dis- 
cussed in subsection III A[ the Hamiltonian for a chain of 
trapped ions coupled by a pair of laser pulses is actually 
Hi(t) [given in Eq. ©], which is related to Hi{t) via the 
phase transformation, ([6]). The solution, \ip(t)) to the 
Schrodinger equation using Hi(t) is therefore related to 
\ip(t)) given above by: 

$(t))=U|V(t)>, (44) 

and so, with the Hamiltonian given in equation ([3]), the 
two preparation steps described in equations (|3 1 [) and 
are simply replaced by 



exp 



| 1 ... 1 ... 0) — ► |0...0)|m). 



N—m 



N 



(45) 



and 

|0...0)|m) 



exp 



iV 



TO|o>. 



(46) 



respectively. Physically, these unwanted phase terms 
arise due to the uneven spacing of the ions along the trap 
axis. They can be corrected by a controlled rotation of 



the state of each ion, j, 
without affecting |0)j or 



so that |1) 



ID 



Alternatively, these con- 



stant phases may be absorbed into the definition of the 
computational basis states. We emphasise that (j) a - ,h are 
known, constant phases, which are unchanged by fluc- 
tuations in the interaction parameters, and that this 
is greatly preferable to the dynamical phases arising in 
many other approaches. 



D. Spatial profile of laser beams and parameter 
fluctuations 



The procedure described above can tolerate signifi- 
cant variations in the experimental parameters. Figure [6] 
makes it clear that the technique is robust against size- 
able fluctuations in Oo and A, and Fig. [7] shows that 
imperfections in the pulse timing have almost no effect 
on the overall fidelity. However, we note that as in pre- 
vious STIRAP-type models, our technique does require 
that the Raman-coupled transition between |0) <-> |1) 
satisfies a condition of two-photon resonance, i.e. 



UJ b + V = CJ 0e - Wi e , 



(47) 



and does not tolerate relative fluctuations in the frequen- 
cies Lo a and ujb- 

We also note that in a practical setting, some slight 
variations in the laser intensity at the position of each 
ion might be expected. In order to quantify the effects of 
such variations, we have simulated the state preparation 
using laser pulses with a Gaussian spatial profile and a 
10% variation in the pulse intensities from the centre to 
the edge of the ion chain. In this case, and using the same 
parameters as in Fig. [5j the overall fidelity remained 
above 98%. 

Finally we remark on another considerable advantage 
of our scheme - namely that the state vector does not ac- 
quire any dynamical or geometric phase; the state \ipo(t)) 
has zero energy in the interaction picture, and the rela- 
tive phases of the laser pulses remain constant during 
the adiabatic passage. The absence of dynamical phases 
represents a significant advantage of our scheme, since 
in other approaches, relatively small unknown imperfec- 
tions in the experimental parameters can lead to sizeable 
unknown phase factors between different components of 
the state vector and this can dramatically reduce the 
overall purity of the state. 



VI. CONCLUSIONS 

We have shown that maximally-entangled Dicke states, 
IH 7 ,^), of arbitrary size may be generated in an ion trap 
using two pairs of extremely simple laser pulses and re- 
maining in a decoherence-free subspace with respect to 
spontaneous emission. Starting with m ions in state |1) 
and all other ions in state |0), a pair of overlapping, 
counter-intuitively ordered pulses - each of which simul- 
taneously addresses these m ions - is used to steer the 
system adiabatically into a phonon number state with 
zero ionic excitations and m phonons. A second pair of 
similar pulses - this time addressing all N ions - is sub- 
sequently used to transfer the system robustly from the 
m-phonon state to the Dicke state, |W^). 

This extremely simple scenario derives from three the- 
oretical findings: (i) The total number of excitations is 
preserved by the Hamiltonian, which is also symmetric 
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with respect to interchange of the ions. For an appropri- 
ately chosen initial state, the dynamics are thus confined 
to the subspace of generalised Dicke states [Eq. (fT2|)]. 
with a given number of excitations; (ii) this symmet- 
ric subspace contains a unique multi-ion dark state, in 
which the upper level |e) is never populated in any ion. 
Therefore, providing that the system remains in this dark 
state, most couplings within the huge symmetric sub- 
space can be neglected and the problem becomes analyt- 
ically tractable; (iii) remarkably, by applying counter- 
intuitively ordered overlapping pulses, this dark state 
can be transformed adiabatically from an initial product 
state into the Dicke state, |W^). 

In addition to an extremely simple experimental im- 
plementation, our proposal also possesses several other 
highly attractive features: (i) the dynamical and geo- 
metric phase acquired during the whole preparation pro- 
cedure are both identically zero; (ii) the proposed tech- 
nique is adiabatic in nature and hence it is robust against 
intensity and frequency imperfections; (iii) there is no de- 
coherence arising from spontaneous emission in the adi- 
abatic limit, regardless of the decay rate from the upper 
level, since the process utilises a dark state. This al- 
lows the use of resonant laser pulses which in turn allows 
shorter pulse durations; (iv) because only two interaction 
steps are required and these can both be performed very 
rapidly, heating effects are greatly reduced when com- 
pared to other schemes. 
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The recurrence relations defining M2 m +i and M 2m are 
therefore: 



M2m+1 — —EtpMim — \ 2 a m M2m-l, 



(A2a) 



M 2m = (A - E v )M 2m _ 1 - \\ m _ 1 M 2m _ 2 , (A2b) 



with 



M x =-E v , (A2c) 
M 2 = — E tp (A — E v ) - , (A2d) 

and elimination of all the cvcn-dimcnsional determinants 
gives: 

M 2m +i = — A^ m M 2m _i — E v (A — E v )M2m-i 



A b.m-1 



(A - E v )M 2m - 3 - \\ m -i 



(A3) 



After recalling that 



A Q:M =~^o -m + n) cxp ( \ 



(A4a) 
(A4b) 
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APPENDIX A: NONZERO EIGENVALUES OF 

Hi(t) 

In addition to the determinant given in Eq. (|37[) . we 
define M 2m as follows: 



M, 



-^6,0 A — E<p \ a ,l 
Xa,l —E v 



-1 A-£ ¥ 



(Al) 



a close inspection of Eq. (|A3[) reveals that the 2m solu- 
tions of Eq. (|36p for which E v ^ are contained in the 
roots of an m -order polynomial, P(z), in the dimen- 
sionless parameter z = E V (A — £ , ¥ ,)/51q. Two adiabatic 
energies, E v , can be derived from each root of P(z). Im- 
portantly for the purposes of subsection IIVB1 the coeffi- 
cients in P(z) are independent of both Qq and A. There- 
fore, the roots of P(z) depend only on the number of 
excitations, m, the number of ions, N, and the time, t. 
We choose the notation z = 7^(t) for the root of P(z) 
corresponding to the energy eigenvalue E±. Therefore, 
the adiabatic energy closest to that of the dark state is: 



E t 



A -a/A 



4«« 



(A5) 
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